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Abstract 

We present an ab initio exploration of the phenomena 
which wiU become important for freestanding struc- 
tures of sihcon as they are reahzed on the nanoscale. 
We find that not only surface but also edge effects 
are important considerations in structures of dimen- 
sions ~3 nm. Specifically, for long nanoscale silicon 
bars, we find two competing low-energy reconstruc- 
tions with a transition from one to the other as the 
cross section of the bar decreases. We predict that 
this size-dependent phase transition has a signature 
in the electronic structure of the bar but little effect 
on elastic properties. 

PACS 68.35.Bs 68.35.Rh 

As our understanding of bulk and surface proper- 
ties of materials matures, the physics of nanoscale 
structures opens new fundamental questions. What 
are the ground state structures of nanoscale collec- 
tions of matter and to what extend can they be pre- 
dicted by simply scaling down bulk and micron-level 
behavior or scaling up the behavior of small clus- 
ters? What new considerations must be taken into 
account? Do nanoscale structures exhibit fundamen- 
tally different electronic or mechanical properties due 
to the large fraction of atoms at surfaces and edges, 
i.e., at the intersection of two surfaces? What are the 
the effects of nanoscale structure on the reconstruc- 
tion of surfaces? 

Clearly, for sufficiently small structures, edges be- 
come important. One key issue is the identification 
of the scale at which this happens and in particular 
whether edges come into play for anything larger than 
a small cluster of atoms. Also, one must determine 
the phenomena by which this importance manifests 
itself. In this letter, we use ab initio calculations to 
show that edge effects indeed become important in 
silicon on length-scales on the order of a few nanome- 
ters, only a factor or two or three times smaller than 



what can be achieved by recent technology [Q. We 
find that the presence of edges has a profound effect 
on the reconstruction on the surface of a structure 
and thereby its electronic structure. Specifically, we 
predict that for long bars along the [001] direction the 
edges drive a surface reconstruction transition from 
the familiar "2x1" family to the "c(2x2)" family[H at 
a cross section of 3 nm x 3 nm. 

Long "bars" (as illustrated in Figure ^) provide the 
ideal laboratory for studying the nature of edges and 
their interaction with surfaces. An isolated edge im- 
plies an infinite system, whereas a bar consists of a 
series of edges bounding a finite area and thus may 
be studied within the supercell framework. In addi- 
tion, such structures are studied experimentally. Us- 
ing lithographic techniques, long bars of silicon can 
be created in the form of suspended bridges between 
bulk silicon supports Q. The heat flow and vibra- 
tional properties of such structures should be unique, 
reflecting quantum conflnement and quantization of 
bulk phononsQ. Furthermore, since the initial re- 
port of bright visible luminescence from "porous 
silicon" Q , there have been many efforts to explain 
this phenomenon based on quantum confinement in 
silicon wires or barsj|, ||, 0, 1). 

In this work, we take on the question of deter- 
mining the ground-state structure of nanometer sized 
bars of silicon as a central issue. Calculations to date, 
where this has not been the central issue, have all 
been done with hydrogen-passivated silicon surfaces, 
placing the atoms at their ideal bulk coordinates |^ or 
simply relaxing to the closest energy minimum with- 
out exploring alternate constructions]^. Hydrogen- 
passivation of silicon surfaces prevents many different 
types of reconstructions, a subject of interest when 
silicon surfaces are exposed to vacuum and which we 
study here. Furthermore, the question of the "round- 
ing" of such bars by the formation of facets along 
their edges has generally been ignored. 

Background information — All of the ab initio elec- 
tronic structure calculations which we report here 
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Figure 1: These two schematic top views of the 
Si(lOO) surface show the dimerization patterns for 
both the 2x1 or p(2x2) (left) and c(2x2) reconstruc- 
tions (right). 




were carried out within the total energy plane wave 
density functional pscudopotential approach[^, us- 
ing the Perdew-Zunger|p^ parametrization of the 
Ceper ly- Alder ijlll exchange-correlation energy and a 
non-local pscudopotential of the Kleiman-Bylander 
forni|jl^ with p and d non-local corrections. In all 
cases, we used a plane wave basis set with a cutoff 
energy of 12 Ry. Electronic minimizations were car- 
ried out using the analytically continued functional 
minimization approach |p^. 

To establish baseline information which we shall 
need later, we computed ab initio energies for var- 
ious reconstructions of the Si(lOO) surface. These 
calculations were carried out in a supercell geome- 
try with slabs of twelve atomic layers containing 48 
silicon atoms and separated by 9A of vacuum. We 
sampled the Brillouin zone using the four k points 
(0,±i,±i). Among the 2x1, p(2x2), and c(2x2) re- 
constructions of the (100) surface (Figure |l|), we find 
in good agreement with previous calculations that 
the p(2x2) is lowest in energy with a binding energy 
of 1.757 eV/dimer, and that the 2x1 reconstruction is 
higher in energy than the p(2x2) by 0.114 eV/dimer. 
Furthermore, we find the c(2x2) reconstruction to be 
higher in energy than the p(2x2) by 0.154 eV/dimer. 
To determine the expected size of the facets along the 
edges of the bars we require the Si(lOO) and Si(llO) 
surface energies. These are known experimentally to 
be 1.36 and 1.43 J m~^ respectively [p^. 

Our ab initio study of edges is carried out on bars 
with a cross section of 2.5 x 2.5 cubic unit-cells in 
the (001) plane. We apply infinite periodic boundary 
conditions in the [001] direction with a periodicity 
of two cubic unit cells. The ball-and-stick diagram 
in Figure |^ depicts the projection of this structure 
on to the (001) plane. Ab initio calculations on the 



Figure 2: Cross sectional view of the silicon bars 
showing the two different Wulff constructions that re- 
sult from using either the experimental (dashed lines) 
or the tight-binding surface energies (dotted Hues). 
The final structure used in this study is formed by 
removing the shaded atoms along the edges. 



bars were performed using the same pscudopotential 
and energy cutoff as used above for the surfaces. For 
the bars, we sampled the Brillouin zone at the two k 
points (0,0, ±i) and provided for a minimum of 6A 
of vacuum between periodic images of the bars. 

Prediction of a finite-size transition — Before de- 
termining the reconstructions along the edges and 
surfaces, we first must determine the overall cross- 
sectional geometry of the bar. To establish this, 
we performed the Wulff construction using the ex- 
perimental surface energies given above. Figure ^ 
shows the resulting shape when {110} facets connect 
{100} surfaces of silicon, scaled to the lateral size of 
our bars. The atomic-scale structure most consistent 
with the Wulff construction in shape and aspect ratio 
appears in Figure ^, where we have removed the four 
columns of shaded atoms along the edges. The pro- 
jection in the (001) plane of the final structure is an 
octagon. With the aforementioned periodicity along 
[001], the final structure contains 114 atoms. With- 
out relaxation, all atoms on the surfaces of the bar 
are two-fold coordinated. 

Having determined the overall geometry, we may 
now turn to the more subtle issue of relaxation and 
possible reconstructions along the surfaces and edges. 
We have found two competing structures: one which 
best satisfies the system in terms of the total num- 
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Figure 3: A schematic view of the bar and the bonds 
formed on the surface of the bar in the c{2x2) recon- 
struction (top) and in the 2x1 reconstruction (bot- 
tom). 

ber of bonds, and the other which best satisfies the 
system in terms of the configuration of the exposed 
surfaces. We find that the system cannot satisfy both 
conditions simultaneously. 

Starting from the unrelaxed configuration, there is 
one unique surface atom with which each edge atom 
may bond in order to become three-fold coordinated. 
This bonding does not change the periodicity along 
the edge, which is the [001] vector of the crystal lat- 
tice. This periodicity, however, is incompatible with 
the periodicity of the p(2x2) low energy state of the 
{100} facets. To maintain maximal bonding, the ar- 
rangement of atoms on the {100} facets must then re- 
vert to the higher energy c(2x2) reconstruction. We 
denote this configuration of the bar as the "c(2x2) 
reconstruction" (See Figure |3|). 

While the c(2x2) reconstruction maximizes the 
number of bonds in the system, the increased energy 
represented by the c(2x2) arrangement on the {100} 
facets will eventually outweigh the benefit of maxi- 
mal bonding along the one dimensional edges for a 
sufficiently large system. When the {100} facets as- 



sume the p(2x2) configuration, the periodicity along 
the [001] direction is doubled. The edge atoms are 
no longer equivalent, and every other atom now can- 
not form a new bond and remains only two-fold co- 
ordinated. Because of bonding restrictions into the 
interior of the structure, the dimer rows along alter- 
nate faces assumes a pattern which leads us to denote 
this configuration as the "2x1 reconstruction." (See 
Figure ||.) 

The ground state of the bar is thus determined by 
the balance between bonding along the edges and the 
surface energy of the facets. We find that for our 
bar the effects from the edges overcome the natural 
tendency of the surfaces. The c(2x2) configuration is 
lower in energy than the 2x1 configuration by 1.94 eV 
per [001] vector along the length of the bar. There 
is therefore a size-dependent transition in this sys- 
tem as we increase the lateral dimension of the bar 
and place relatively more atoms on the {100} facets. 
Our calculations place the crossover point at approx- 
imately 5 dimer pairs on each (100) surface per unit 
cell along [001], corresponding to a bar with cross 
sectional dimensions of approximately 3.0 nm x 3.0 
nm, far larger than the scale of an atomic cluster. 
We therefore predict an important edge-driven tran- 
sition at a scale only two or three smaller than what 
has been achieved experimentally to date[|l|. 

More generally, we can see that the compatibility of 
the competing surface reconstructions with the trans- 
lational symmetry of the edges plays a pivotal role in 
determining the ground state of nanoscale structures. 
In our specific case, it is the principal physical mech- 
anism giving rise to the size-dependent transition for 
our bars. 

Implications of the transition — Comparing the 
electronic structures of the two reconstructions (see 
Figure ^), we find that the 2x1 configuration has a 
gap of 0.35 eV across the Fermi level. Furthermore, 
the topmost filled states consist of four nearly de- 
generate states which are localized on the four edges 
of the bars, and this cluster is separated from states 
below it in energy by a gap of 0.30 eV. This nearly 
symmetric placement with sizeable gaps as well as 
the spatial localization of the edge states leads us to 
conclude that the 2x1 bar is insulating. 

On the other hand, the electronic structure of the 
c(2x2) configuration is more subtle. The states in 
the vicinity of the Fermi level are localized on the 
surfaces of the bars, and the gap across the Fermi 
level is only 0.09 eV. We believe the c(2x2) configu- 
ration is a small-gap semiconductor or even perhaps 
metallic. Thus the nature of the low energy electronic 
states and excitations differ between the two bars and 
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Figure 4: Eigen-energies of the electronic states for 
the two reconstructions of the bars in the vicinity of 
their Fermi levels at the k points used in the calcula- 
tions: (a) is the 2x1 configuration, (b) is the c(2x2) 
configuration. Filled states are denoted by filled cir- 
cles and empty states by empty circles. The zero of 
energy is arbitrary. 

should manifest themselves in physical measurements 
such as electrical conductivity or optical spectra. 

Next, considering possible differences in mechan- 
ical properties, we carried out molecular dynamics 
simulations using a tight-binding model to compute 
transverse acoustic phonon frequencies for the bars. 
We believe that the tight-binding model will give us 
a qualitative view of the differences which may ex- 
ist between the two bars. If any large differences are 
found, we should examine this issue in more detail 
using the more demanding ab initio techniques. 

We used the semi-empirical tight-binding model 
of Sawad a[|l6| with the modification proposed by 
Kohyama|l7|. This model provides a good qualita- 
tive description of bulk, dimer, and surface energetics 
of silicon. In using the Sawada model, we only keep 
Hamiltonian matrix elements and repulsive terms be- 
tween atoms closer than r„„ = 6 A. 

In order to calculate the band-structure energy, 
we used the fully parallelizable 0{N) technique of 
Goedecker and Colombo |l8|. By checking the con- 
vergence of the total energy to its ground-state value 
(determined by exact diagonalization), we found it 
necessary to use the following set of parameters to 
ensure convergence of energies to within an accuracy 
of 10^"* eV/atom: in the nomenclature of we 
have fcsT = 0.125 eV, Upi = 300, and rioc = 15.0 A. 

Table summarizes the results of the Sawada 
model for the various key physical values in this 
study. Although the Sawada model correctly predicts 
asymmetric dimers as the ground-state of the Si(lOO) 
surface, it is not sensitive to the delicate rocking of 
dimers that differentiates the 2x1 and p(2x2) recon- 
structions and hence finds the 2x1 reconstruction to 
be lower in energy than the p(2x2). However, in our 



study, the relevant energy difference is that between 
the c(2x2) and the lowest-energy reconstruction of 
the surface, and this quantity is reproduced rather 
well. The Sawada model also does well in predict- 
ing surface energies: Figure ^ shows the result for 
the Wulff construction when we use the tight-binding 
surface energies, and the resulting geometry is very 
similar to the experimentally derived one. Thus we 
believe that the Sawada model provides a good semi- 
quantitative description for the physics of our bars. 

Using the Sawada model, we computed transverse 
acoustic phonon frequencies for the k = (7ra/4)i 
mode (along A) for both reconstructions of our bars, 
the longest allowed wavelength along the length of 
the bar consistent with the periodic boundary condi- 
tions. Wc ran {N,V,E) molecular dynamics simula- 
tions using the Verlet algorithm with a time-step of 
2.4 fs for 900 time steps. Phonon frequencies were 
identified as peaks in the frequency-domain power- 
spectrum of the velocity autocorrelation function as 
estimated by auto-regressive fits. We found frequen- 
cies of 1.98 ± 0.02 THz and 1.93 ± 0.02 THz for the 
c(2x2) and 2x1 configurations respectively. Thus we 
see that edge effects do not have a significant effect 
on the long wavelength vibrations of the bars. 

Conclusions — We have performed an ab initio 
study of the energetics of long bars of silicon in vac- 
uum. We found useful the atomic-scale version of 
the Wulff construction as a first step in determining 
nanoscale structure. Next, we found that one can- 
not ignore the interplay between the edges and sur- 
faces in silicon structures with dimensions of a few 
nanometers, where the compatibility of the surface 
reconstructions with the symmetry of the edges plays 
an important role. In particular, the ground-state of 
our bars changes from one surface reconstruction to 
another, signalling a cross section-dependent phase 
transition, with significant influence on the electronic 
structure of the system. Finally, we find that even on 
the scale of a few nanometers, the surface-edge inter- 
play has little effect on mechanical properties. 
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Expt/Ab initio 


TB 


Bulk properties 


Binding energy (eV/atom) 


-4.63" 


-4.79 


Bulk modulus (Mbar) 


0.975'^ 


0.906 


Phonon frequencies 


r (THz) 


15.5" 


18.1 


k = (7r/4a)z LA (THz) 


3.9'' 


3.8 


k = (7r/4a)i TA (THz) 


2.4'' 


3.0 


Si(lOO) reconstructions 


lowest energy (eV/dimer) 


p(2x2): 1.757 


2x1: 2.05 


c(2x2) (eV/dimer) 


1.603 


1.93 


Surface energies 


Si(lOO) (J m-^) 


1.36= 


1.44 


Si(llO) (J m-^) 


1.43^= 


1.77 


Differences between zizag and c(2x2) bars 


Energy difference (eV/a) 


1.94 


1.57 


Cross-over (width in nm) 


3.0 


3.0 



Table 1: Comparison of the Sawada model (TB) with experimental and ab initio results: (a) is reference 
| p9[ , (b) is reference and (c) is reference |l^. The energy difference between the 2x1 and c(2x2) 

reconstructions of the silicon bars are given in eV per unit cell along [001], and the cross-over refers to the 
approximate width of a bar when the two reconstructions have the same energy (see text). 
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